Optical bistability at low light level due to collective atomic recoil 
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We demonstrate optical nonlinearities due to the interaction of weak optical fields with the col- 
lective motion of a strongly dispersive ultracold gas. The combination of a recoil-induced resonance 
(RIR) in the high gain regime and optical waveguiding within the dispersive medium enables us to 
achieve a collective atomic cooperativity of 275 ± 50 even in the absence of a cavity. As a result, we 
observe optical bistability at input powers as low as 20 pW. The present scheme allows for dynamic 
optical control of the dispersive properties of the ultracold gas using very weak pulses of light. The 
experimental observations are in good agreement with a theoretical model. 



Motivated by potential applications to quantum infor- 
mation science jl|, 0, 0, there have been intense ex- 
perimental efforts to realize strong nonlinear interactions 
between dilute atomic ensembles and weak optical fields. 
Methods to achieve such quantum nonlinear couplings 
in dissipation-free media have relied mainly on two ap- 
proaches. First, the interaction time and the coupling 
between the atoms and the photons can be enhanced by 
placing the atoms within a high finesse cavity [5j , an ap- 
proach that comes at the expense of considerable exper- 
imental complexity and low bandwidth. Alternatively, 
near-resonant light propagating through optically dense 
media can also result in strong nonlinear interactions. In 
order to limit the dominant linear absorption, this lat- 
ter approach requires the use of a coherent multiphoton 
process such as EIT An added benefit of such a 
multiphoton process is the enhancement of the interac- 
tion time due to slow group velocities. However, since 
EIT relies on quantum interference between the internal 
states of an atom, it is fairly sensitive to inhomogenous 
optical and magnetic fields. 

In this Letter, we demonstrate nonlinear optical effects 
due to the interaction between weak pulses of light and 
the collective motion of an ultracold slow-light medium. 
Using the motional degrees of freedom to create a highly 
dispersive gas alleviates the sensitivity to external fields. 
As shown in recent studies 0, HI , the implementation of 
a recoil-induced resonance (RIR) in an optically dense 
anisotropic gas allows for strong atom-light interaction 
due to the combination of slow group velocities and trans- 
verse confinement of the optical fields within the atomic 
medium. Due to this strong coupling, we observe optical 
bistability at input powers as low as 20 pW, an upper 
bound limited mainly by the photodetector efficiency. 

The motion of delocalized atoms under the influence 
of two or more light fields can mediate the conversion of 
atomic kinetic energy into radiation [9(. These processes, 
termed 'recoil-induced resonances', can be described in 
terms of stimulated Raman transitions between different 
momentum classes of the atomic ensemble (lo| . For an 
optically thin medium, the atom-light interaction has lit- 
tle effect on the momentum distribution of the gas and 
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FIG. 1: (a) Schematic indicating the pump and probe beams 
used for the RIR. (b) The energy levels relevent to the RIR. 
(c) Absorption spectrum of a probe beam around the RIR for 
longitudinal optical densities of ~ 10 (i) and ~ 40 (ii). 



a perturbative analysis reveals that the probe field ex- 
periences weak absorption (gain) for positive (negative) 
detuning relative to the pump field. 

In contrast, as the optical density of the atomic gas 
increases, the strong atom-light coupling leads to large 
optical gain Q and significant effects on the momentum 
distribution of the atomic gas. The strong amplification 
of the probe field and its subsequent back-action on the 
collective motional states of the gain medium results in a 
nonlinear optical response even for weak incident beams. 
In this work, this back-action is evidenced by the obser- 
vation of optical bistability in the probe transmission. 

The experimental scheme to create a strongly disper- 
sive ultracold gas using the RIR is described in previous 
work Q. About 5 x 10 s 87 Rb atoms are confined in a 
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FIG. 2: (a) Absorption bistability due to the interaction be- 
tween the probe and the collective motion of the atomic gas. 
Depending on the sign of the detuning chirp (dS/dt), the 
transmission indicates a shift in both the resonance and the 
peak amplification (g+,g~). (b) Peak probe amplification in- 
dicates a bistable response down to the detection limit of the 
input probe power. 



highly anisotropic magneto-optic trap at typical temper- 
atures of ~ 20 fiK. In this trap, the atom cloud assumes 
the shape of a cylinder with approximate radial (longi- 
tudinal) extent of 200 /xm (3 cm). The pump and probe 
beams for the RIR share the same linear polarization 
and are directed along the long axis of this cylinder to 
take advantage of the large optical depth (OD) along this 
axis. The Rabi frequency of the pump beam was typi- 
cally Qi/r = 1.5 where T is the natural linewidth of this 
transition. At low OD, the transmission spectrum of the 
RIR exhibits a characteristic dispersion-shaped spectrum 
with gain (absorption) for detuning d < (6 > 0). At 
higher OD, the probe is almost completely extinguished 
on the absorption side of the resonance leaving the gain 
peak as the only distinguishable feature (Fig. 1). 

For an atomic gas with high OD, Fig. 2(a) shows 
the probe transmission as the pump-probe detuning is 
scanned across the RIR. Depending on the sign of this 
detuning chirp, there a shift in the resonance as well 
as in the maximum gain. This hysteretic nature of the 
transmission diminishes as the OD is lowered to less than 
~ 10. At larger cooperativity, the atomic ensemble and 
the light fields form a strongly coupled system and the 



probe interacts with an atomic ensemble whose motional 
coherences and momentum distribution is the cumulative 
result of prior interactions with the light fields. As seen 
in Fig. 2(b), the bistable response (g~/g+ > 1) persists 
down to the detection limit (~ 20 pW) of the input probe 
power. 

In order to understand this behavior, we first note that 
the amplification of the probe is accompanied by the 
transfer of atoms from a momentum p to a momentum 
p + 2hk (Fig. 1(b)). Thus, as the detuning is scanned 
across the gain side of the RIR, atoms at various mo- 
menta are brought into resonance with the light fields 
and transferred to higher momentum states. Crucially, 
in the case of a negative chirp (dS/dt < 0), these trans- 
ferred atoms are brought closer to resonance with the 
light fields as the detuning is scanned. Thus, in this 
process, atoms are progressively swept up the momen- 
tum ladder due to the time- varying detuning, resulting in 
both an enhanced amplification and a shift in the location 
of the RIR. In contrast, for a positive chirp (dS/dt > 0), 
the atoms are transferred to states that are farther from 
resonance and the probe transmission resembles that ob- 
tained for a static thermal distribution. 

To quantitatively explain these observations, we use a 
theoretical model that describes two classical light fields 
that are coupled to the motional degrees of freedom of an 
elongated ensemble of two-level atoms. To mimic the ex- 
periment, the atoms are assumed to be tightly confined 
in the radial dimension. Accordingly, only the atomic 
momentum along the lon g ax is is relevent and the Hamil- 
tonian can be written as 
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where c g (k)(c e (k)) are the annihilation operators of 
ground (excited) state atoms with momentum hk, ujq is 
the transition frequency of the two-level atoms, 31(172) is 
the atom- light coupling coefficient and 01(02) is the nor- 
malized electric field of the pump (probe) beams. In the 
far-detuned limit, the excited states can be adiabatically 
eliminated and the equation for the coherences and pop- 
ulations of the different momentum classes in the ground 
state arc 
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P(p,p') = 4iw r (p' 2 ~p 2 )p(p,p') 

a 2 e- l5t (p(p+l,p>)-p(p,p> -1)) 

(-p(p,p'+l)+p(p-l,p')) (2) 

Uk') 



A z 

where p(p,p') = (c g (k')c g (k)), ui r is the recoil frequency 
and S = LU2 — uj\ is the detuning. Also, using the slowly- 
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varying envelope and single- mode approximations, the 
evolution of the probe amplitude can be written as 



dt 



a 2 = iN 



ai e t5t J2p(P~hp)~^a 2 -a m ). (3) 



Retaining only first-order coherence terms between 
momentum classes, the above equations can be writ- 
ten as a set of coupled equations for the population 
flp = p(p,p), the first-order coherence rj p = p(p+l,p)e 
and the probe amplitude a 2 . 
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where 7 p0 p(7coh) are the population (coherence) relax- 
ation rates respectively and p = gig 2 ai/A. The ther- 
mal population distribution Hth,p is given by Maxwell- 
Boltzman distribution. The decay rate of photons is 
approximated by the free-space rate n = c/L with L 
the longitudinal extent of the atomic gas. These cou- 
pled equations describe the rich dynamics that ensue as 
a consequence of the collective atom-light interaction and 
a time-dependent pump-probe detuning. 

Numerical simulations of the bistability based on this 
model show excellent agreement with the experimental 
results over a wide range of parameters of pump detuning 
and chirp rate (Fig. 3). In these simulations, parameters 
such as the pump detuning, OD of the atomic gas and 
scan rate of the two-photon detuning were held fixed at 
the experimental values. 

The observation of a bistable transmission requires 
that the momentum coherences established as a result 
of the atom-light interaction persist as the detuning is 
scanned across the RIR. In practice, off-resonant light 
scattering causes the thermalization of the momentum 
distribution thereby suppressing the bistability. This 
competing process leads to a limiting rate for the scan 
rate (dd/dt) below which the atomic momentum distri- 
bution is quasi-static and the transmission becomes non- 
hysteretic (Fig. 4(a)). A bistable transmission was ob- 
served for scan rates as low as 500 kHz/ms, correspond- 
ing to moving across the RIR in 100 /is, much longer than 
the mean photon scattering time of < 1 /is. This obser- 
vation of long-lived momentum coherences is consistent 
with previous studies of the RIR 0, HH • 

An independent measure of the influence of off- 
resonant scattering on the probe transmission was ob- 
tained by determining the thermalization time (7p~ p) 
of the atomic momentum distribution. For this, the 
probe was switched on for 100 ms at an intensity of 
0.1 mW/cm 2 and at a two-photon detuning that cor- 
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FIG. 3: Numerical simulations based on the theoretical model 
match the experimental results over a wide range of parame- 
ters. Panels on the left show the numerical simulations for the 
gain coefficient given by exp(— 2Re[a]L) where a is the ab- 
sorption coefficient and panels on the right show the observed 
experimental gain of the probe transmission across the RIR, 
both versus two-photon detuning. (a), (b) and (c) correspond 
to transmission spectra obtained by chirping the two-photon 
detuning at scan rates of 0.1, 0.5 and 2.5 MHz/ms, respec- 
tively. 



responded to the peak gain of the RIR. The probe in- 
tensity was then reduced to ~ 10~ 3 mW/cm 2 within a 
few microseconds. Following this reduction in intensity, 
the probe transmission indicated a gain that was initially 
very small but gradually relaxed to a higher equilibrium 
value. This relaxation was interpreted as being due to 
the thermalization of the momentum distribution to re- 
populate those states that were depleted due to the ini- 
tially intense probe. Under the assumption that the fi- 
nal probe intensity was too small to significantly modify 
the atomic distribution, the time scale over which the 
probe transmission reaches a steady state should reflect 
the thermalization time. Consistent with the observa- 
tion of bistability at relatively slow scan rates, we mea- 
sure thermalization times on the order of a few hundred 
microseconds (Fig. 4(b)). 

The numerical simulations used the measured value of 
Ipop (Fig. 4(b)) while allowing the value of j co h to vary. 
In order to obtain the best match with the experimental 
results, we typically found that ^ C oh was required to be 
greater than 7p 0p . However, both quantities were much 
smaller than the photon scattering rate. This sup pression 
is in agreement with earlier observations 0, [II, HI, 14 1 ■ 
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A possible explanation lies in the fact that in the pres- 
ence of an optical potential (U oc 1/A), in the limit of 
large detuning, both 7^ and 7 pop can be suppressed 
due to Lamb-Dicke confinement. The estimated scaling 
of the suppressed decay rate with the detuning should be 
j e ff oc ^scatt/VU oc A~ 3/2 where j sca tt is the photon 
scattering rate. The experimental data (Fig. 4(b)) also 
suggests a similar scaling(1.57 ± 0.09). 

The observation of optical bistability at low input pow- 
ers suggests prospects of all-optical control using weak 
pulses of light. For instance, this scheme lends itself 
rather easily to a low light level all-optical switch wherein 
the detuning of a few-photon probe is controlled in or- 
der to switch a more intense output beam. The largest 
scan rates, dS/dt = 30 MHz/ms at which bistability 
was observed and typical values of the momentum co- 
herence time y~J~ h ~ 100 fj,s together yield an estimate of 
the switching time t = jcoh/ (dS/dt) — 0.3 fj,s. This is 
commensurate with previous measurements of all-optical 
switching using a MR [7]. Combining this with the low- 
est probe powers (P ~ 20 pW) and typical beam waists 
(100 /im) used in this work, we obtain a typical photon 
number tP/ (hc/X) ~ 25 and a remarkably low switching 
energy density [3, of 7 x 10~ 5 photons/(A 2 /27r) to 
operate this all-optical switch. 

In conclusion, we demonstrate optical nonlinearities 
due to a coherent interaction between weak light fields 
and the collective motion of a strongly dispersive atomic 
gas. Since the atomic momentum is relatively insensitive 
to external magnetic or electric fields, such systems may 
be promising candidates for applications in low-light level 
nonlinear optics. 
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FIG. 4: (a) The ratio of peak gain (g-/g+) indicates a limiting 
scan rate dS/dt ~ 0.5 MHz/ms below which decoherence and 
thermalization of the momentum distribution suppress opti- 
cal bistability. This data was obtained at a pump detuning 
A ~ — 4T. (b) Thermalization time (7^) vs the pump de- 
tuning. A fit to the data (dashed line) indicates a power-law 
dependence 7j7 p oc A" with a — 1.57 ± 0.09. 



